Cosmological Singularities, Einstein Billiards 
and Lorentzian Kac-Moody Algebras 

Thibault Damour 

Institut des Hautes Etudes Scientifiques, 35 route de Chartres, 
91440 Bures-sur-Yvette, France 



Abstract 

The structure of the general, inhomogeneous solution of (bosonic) 
Einstein-matter systems in the vicinity of a cosmological singular- 
ity is considered. We review the proof (based on ideas of Belinskii- 
Khalatnikov-Lifshitz and technically simplified by the use of the Arnowitt- 
Deser-Misner Hamiltonian formalism) that the asymptotic behaviour, 
as one approaches the singularity, of the general solution is describ- 
able, at each (generic) spatial point, as a billiard motion in an aux- 
iliary Lorentzian space. For certain Einstein-matter systems, notably 
for pure Einstein gravity in any spacetime dimension D and for the 
particular Einstein-matter systems arising in String theory, the bil- 
liard tables describing asymptotic cosmological behaviour are found 
to be identical to the Weyl chambers of some Lorentzian Kac-Moody 
algebras. In the case of the bosonic sector of supergravity in 11 di- 
mensional spacetime the underlying Lorentzian algebra is that of the 
hyperbolic Kac-Moody group E\q, and there exists some evidence of 
a correspondence between the general solution of the Einstein-three- 
form system and a null geodesic in the infinite dimensional coset space 
E\q/K{E\q), where K{E\q) is the maximal compact subgroup of E\q. 



1 Introduction and general overview 

A series of works [24, 25, 26, 27, 23, 28, 29] has uncovered a remarkable con- 
nection between the asymptotic behaviour, near a cosmological singularity, 
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of certain Einstein-matter systems and billiard motions in the Weyl chambers 
of some corresponding Lorentzian Kac-Moody algebras. This simultaneous 
appearance of billiards (with chaotic properties in important physical cases) 
and of an underlying symmetry structure (infinite-dimensional Lie algebra) is 
an interesting fact, which deserves to be studied in depth. Before explaining 
in detail the techniques (including the Arnowitt-Deser-Misner Hamiltonian 
formalism [1]) that have been used to uncover this fact, we shall start by 
reviewing previous related works, and by stating the main results of this 
billiard/symmetry connection. 

The simplest example of this connection concerns the pure Einstein sys- 
tem, i.e. Ricci-flat spacetimes. Long ago, Belinskii, Khalatnikov and Lif- 
shitz (BKL) [8, 9, 10] gave a description of the asymptotic behaviour, near 
a spacelike singularity, of the general solution of Ricci tensor = in (3 + 1)- 
dimensional spacetime in terms of a continuous collection of second-order, 
non linear ordinary differential equations (with respect to the time vari- 
able). As argued by these authors, near the singularity the spatial points 
x l , % — 1,2, 3, essentially decouple, and this decoupling allows one to approx- 
imate a system of 'partial differential equations (PDE's) in 4 variables (t, x % ) 
by a 3-dimensional family, parametrized by (x l ) G M 3 , of ordinary differen- 
tial equations (ODE's) with respect to the time variable t. The coefficients 
entering the nonlinear terms of these ODE's depend on the spatial point 
x % but are the same, at each given x\ as those that arise in some spatially 
homogeneous models. In the presently considered 4-dimensional "vacuum" 
(i.e. Ricci-flat) case, the spatially homogeneous models that capture the 
behaviour of the general solution are of the Bianchi type IX or VIII (with 
homogeneity groups SU(2) or SL(2,R), respectively). The asymptotic evo- 
lution of the metric was then found to be describable in terms of a chaotic 
[60, 18] sequence of generalized Kasner (i.e. power-law) solutions, exhibit- 
ing "oscillations" of the scale factors along independent spatial directions 
[8, 9, 10], or, equivalently, as a billiard motion [19, 64] on the Lobachevskii 
plane. 

This picture, found by BKL in the 4-dimensional pure Einstein case, was 
then generalized in various directions: namely by adding more spatial di- 
mensions, or by adding matter fields. Remarkably, it was found that there 
exists a critical spacetime dimension for the asymptotic behaviour of the gen- 
eral Ricci-flat solution [34]. When D < 10 the solution exhibits a BKL-type 
never-ending oscillatory behaviour with strong chaotic properties, while when 
D > 11 the general solution ceases to exhibit chaotic features, but is instead 
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asymptotically characterized by a monotonic Kasner-like solution [34, 30]. 
The addition of matter fields was also found to feature a similar subcriti- 
cal/overcritical classification [7, 24, 30]. Let us consider, in some given space- 
time dimension D, the Einstein-dilaton-p-form system, i.e. let us add a mass- 
less scalar 0, and one or several p-form fields A = A Ml ... Mp dx^ 1 A ... A dx^ p /p\. 
It was found that the general solution, near a cosmological singularity, of 
this system is monotonic and Kasner-like if the dilaton couplings X p (defined 
below) of the p-forms belong to some (dimension dependent) open neighbour- 
hood of zero [30] , while it is chaotic and BKL-like if the dilaton couplings be- 
long to the (closed) complement of the latter neighbourhood. In the absence 
of any p-forms, the Einstein-dilaton system is found to be asymptotically 
monotonic and Kasner-like [7, 2]. 

A convenient tool for describing the qualitative behaviour of the general 
solution near a spacelike singularity is to relate it to billiard motion in an 
auxiliary Lorentzian space, or, after projection on the unit hyperboloid, to 
billiard motion on Lobachevskii space. This billiard picture naturally follows 
from the Hamiltonian approach to cosmological behaviour (initiated long ago 
[63] in the context of the Bianchi IX models; see [49] for a recent review), 
and was first obtained in the 4-dimensional case [19, 64] and then extended 
to higher spacetime dimensions with p-forms and dilaton [54, 56, 45, 46, 26]. 
Recent work has improved the derivation of the billiard picture by using 
the Iwasawa decomposition of the spatial metric and by highlighting the 
general mechanism by which all the "off-diagonal" degrees of freedom (i.e. 
all the variables except for the scale factors, the dilaton and their conjugate 
momenta) get "asymptotically frozen" [29]. 

The billiard dynamics describing the asymptotic cosmological behaviour 
of D-dimensional Einstein gravity coupled to n dilatons <f>i, . . . ,<f> n (having 
minimal kinetic terms — Ej \d<f)j\ 2 ), and to an arbitrary menu of p-forms (with 
kinetic terms ~ — e^ Xp '^ \dA p \ 2 , where \ p is a linear form) consists of a "billiard 
ball" moving along a future-directed null geodesic (i.e. a "light ray") of a 
(D — 1 + n)-dimensional Lorentzian space, except when it undergoes specular 
reflections on a set of hyperplanar walls (or "mirrors" ) passing through the 
origin of Lorentzian space. One can refer to this billiard as being an "Einstein 
billiard" (or a "cosmological billiard"). 

As we shall see below in detail the degrees of freedom of the billiard ball 
correspond, for D — 1 of them, to the logarithms, — 2f3 a ) of the scale factors, 
i.e. the diagonal components of the (D — 1) x (D — 1) spatial metric in 
an Iwasawa decomposition, and, for the remaining ones, to the n dilatons 
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(f)j. The walls come from the asymptotic elimination (in the Arnowitt-Deser- 
Misner Hamiltonian [1]) of the remaining "off-diagonal" degrees of freedom: 
off-diagonal components of the metric, p-form fields, and their conjugate 
momenta. It is often convenient to project the Lorentzian billiard motion 
(by a linear projection centered on the origin) onto the future half of the 
unit hyperboloid G ^ ^ f3 u = — 1, where {(3^) = {(3 a , tfij) denotes the position 
in, and Gr M „ the metric of (D — 1 + n)-dimensional Lorentzian space. [The 
signature of is — h + . . . +.] This projection leads to a billiard motion in 
a domain of hyperbolic space Hc-2+n, see Fig. 1. This domain (the "billiard 
table" x ) is the intersection of a finite number of half hyperbolic spaces, say the 
positive sides of the hyperbolic hyperplanar walls defined by intersecting the 
original Lorentzian hyperplanar walls with the unit hyperboloid G^ v /3 U = 
-1. 

The distinction mentioned above between a chaotic, BKL-like behaviour 
and a monotonic, Kasner-like one then corresponds to the distinction between 
a finite-volume billiard table (in H D _ 2 +n) and an infinite-volume one. One 
could further distinguish: (i) the overcritical case {compact billiard table), 
(ii) the critical case [finite-volume, but non-compact billiard table) and (iii) 
the subcritical case {infinite-volume billiard table). Both cases (i) and (ii) 
lead to never-ending, chaotic oscillations of the scale factors (and of the 
dilatons), while case (iii) leads, after a finite number of reflections off the 
walls {i.e. after a finite number of oscillations of the scale factors and of the 
dilatons), to an asymptotically free motion of /9 M (null geodesic), correspon- 
ding to an asymptotically monotonic, Kasner-like (power-law) behaviour of 
the spacetime metric near the singularity, see Fig. 1. In actual physical 
models, only the critical and subcritical cases are found. 

Having introduced the concept of (Lorentzian or, after projection, hyper- 
bolic) billiard table associated to general classes of Einstein-matter systems, 
we can now describe the connection uncovered in [24, 25, 26, 27, 23, 28, 29] 
between certain specific Einstein-matter systems and Lorentzian Kac- Moody 
algebras. In the leading asymptotic approximation to the behaviour near the 
cosmological singularity, this connection is simply that the (Lorentzian) bil- 
liard table describing this behaviour can be identified with the (Lorentzian) 
Weyl chamber of some corresponding (model-dependent) Lorentzian Kac- 

1 We use the term "billiard table" to refer to the domain within which the billiard 
motion is confined. Note, however, that this "table" is, in general, of dimension higher 
than two, and is either considered in Lorentzian space or, after projection, in hyperbolic 
space. 
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Figure 1: Sketch of billiard tables describing the asymptotic cosmo logical 
behaviour of Einstein-matter systems. The upper panels are drawn in the 
(D— l+n)-dimensional auxiliary Lorentzian space spanned by = (j3 a , <ft). 
[D — 4 and n = in the case drawn here.] The billiard tables are represented 
as "wedges" in [(2 + 1) -dimensional] /3-space, bounded by hyperplanar walls 
wa{/3) = on which the (unrepresented) billiard ball undergoes specular 
reflections. The upper left panel is a (critical) "chaotic" billiard table (con- 
tained within the /3-space future light cone), while the upper right one is a 
(subcritical) "non-chaotic" one (extending beyond the light cone). The lower 
panels represent the corresponding billiard tables (and billiard motions) af- 
ter projection onto hyperbolic space [if 2 in the case drawn here]. The latter 
projection is defined in the text by central projection onto 7-space (i.e. the 
unit hyperboloid G^^^f = —1, see upper panels), and is represented in 
the lower panels by its image in the Poincare ball [disc]. "Chaotic" billiard 
tables have finite volume in hyperbolic space, while "non-chaotic" ones have 
infinite volume. 
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Moody algebra. We recall that the Weyl chamber of a Lie algebra with simple 
roots 2 cej, % — 1, . . . , r (r being the rank of the algebra) is the domain of the 
Cartan subalgebra (parametrized by (3 G W) where (ctj, f3) > for all i's. For 
this connection to be possible many conditions must be met. In particular: 
(i) the billiard table must be a Coxeter polyhedron, i.e. the dihedral angles 
between adjacent walls must be integer submultiples of n (i.e., of the form 
7T jk where k is an integer > 2), and (ii) the billiard table must be a simplex, 
i. e. have exactly D — 1 + n faces. It is remarkable that this seemingly very 
special case of a "Kac-Moody billiard" is found to occur in many physically 
interesting Einstein-matter systems. For instance, pure Einstein gravity in 
Z}- dimensional spacetime corresponds to the Lorentzian Kac-Moody algebra 
AEu-i [27]. The latter algebra is the canonical Lorentzian extension [52] 
of the ordinary (finite-dimensional) Lie algebra Ads (= si 0-2)- Another 
interesting connection between qualitative PDE behaviour and Kac-Moody 
theoretic concepts is that the transition between "critical" chaotic, BKL- 
like behaviour in "low dimensions" (D < 10) and "subcritical" monotonic, 
Kasner-like behaviour in "high dimensions" (D > 11) is found to be in strict 
correspondence with a transition between an hyperbolic Kac-Moody algebra 
(AEd-i for D < 10) and a non hyperbolic one (AEd-i for D > 11). We 
recall here that V. Kac [52] defines a Kac-Moody algebra to be hyperbolic 
by the condition that any subdiagram obtained by removing a node from its 
Dynkin diagram be either of finite or affme type. Hyperbolic Kac-Moody 
algebras are necessarily Lorentzian (i.e. the symmetrizable Cartan matrix is 
of Lorentzian signature), but the reverse is not true in general. 

Another connection between physically interesting Einstein-matter sys- 
tems and Kac-Moody algebras concerns the low-energy bosonic effective ac- 
tions arising in String and M-theory. Bosonic String theory in any space- 
time dimension D is related to the Lorentzian Kac-Moody algebra DEd 
[26, 23]. The latter algebra is the canonical Lorentzian extension of the 
finite-dimensional algebra -Dd-2- The various Superstring theories (in the 
critical dimension D = 10) and M-theory have been found [26] to be related 
either to E w (when there are two sypersymmetries in D = 10, i.e. for type 
IIA, type IIB and M-theory) or to BE 10 (when there is only one supersym- 
metry in D = 10, i.e. for type I and the two heterotic theories), see Fig. 2. 

2 The roots a are linear forms on a Cartan subalgebra t) that appear as eigenvalues of 
the adjoint action of t) on the rest of the Lie algebra: [h, e a ] — a(h)e a for all h £ f) . In 
the simple case of the Lie algebra of SU(2), \) — { J z } is one-dimensional and the "raising" 
or "lowering" operators e a — J x ± iJ y correspond to the roots a(J z ) = ±1. 
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Figure 2: Coxeter-Dynkin diagrams encoding the geometry of the billiard 
tables describing the asymptotic cosmological behaviour of three blocks of 
string theories: B2 = {M-theory, type IIA and type IIB superstring theories}, 
£>i = {type I and the two heterotic superstring theories}, and £>o = {closed 
bosonic string theory in D = 10}. Each node of the diagrams represents 
a dominant wall of the corresponding cosmological billiard. Each Coxeter 
diagram of a billiard table corresponds to the Dynkin diagram of a correspon- 
ding (hyperbolic) Kac- Moody algebra: E 10 , BE W and DE 10 , respectively. In 
other words, the cosmological billiard tables can be identified with the Weyl 
chambers of the corresponding Lorentzian Kac-Moody algebras. 

See Ref. [40] for a discussion of the billiards associated to all types of pure 
super gravities in D = 4. See Ref. [23] for the construction of Einstein-matter 
systems related (in the above "billiard" sense) to the canonical Lorentzian 
extensions of all the finite-dimensional Lie algebras (A n , B n , C n , D n , G2, F4, 
E 6 , E 7 and E s ). See also Ref. [32] for the identification of all the hyperbolic 
Kac-Moody algebras whose Weyl chambers are Einstein billiards. 

Ref. [28] has studied in more detail the correspondence between the spe- 
cific Einstein-three-form system (including a Chern-Simons term) describ- 
ing the bosonic sector of 11-dimensional supergravity (also known as the 
"low-energy limit of M-theory") and the hyperbolic Kac-Moody group E w . 
Ref. [28] introduced a formal expansion of the field equations in terms of 
positive roots, i.e. combinations a = T li n l ai of simple roots of E w , a*, 
i = 1, . . . , 10, where the n*'s are integers > 0. It is then useful to order this 
expansion according to the height of the positive root a = ^n 1 cti, defined 
as ht(a) = Ejn\ The correspondence discussed above between the leading 
asymptotic evolution near a cosmological singularity (described by a billiard) 
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and Weyl chambers of Kac-Moody algebras involves only the terms in the 
field equation whose height is ht(a) < 1. By contrast, Ref. [28] could show, 
by explicit calculations, that there exists a way to define, at each spatial 
point x, a correspondence between the field variables g^it, x), A^ u \(t,x) 
(and their gradients), and a (finite) subset of the parameters defining an el- 
ement of the (infinite-dimensional) coset space E w / K(E W ) [where K(E 10 ) 
denotes the maximal compact subgroup of E w ], such that the (PDE) field 
equations of supergravity get mapped onto the (ODE) equations describing 
a null geodesic in E w / K(E W ) up to terms of height 30. 

A complementary check of the correspondence between 11-dimensional 
supergravity and the Ei /K(Ei ) (x-model has been obtained in [31]. This 
result was further extended to the correspondence between the E w / K(E W ) 
cr-model and, both, massive 10-dimensional IIA supergravity [57] and ten- 
dimensional IIB supergravity [58]. 

These tantalizing results suggest that the infinite-dimensional hyperbolic 
Kac-Moody group E 10 may be a "hidden symmetry" of supergravity, in the 
sense of mapping solutions onto other solutions (the appearance of E w as 
possible symmetry group of supergravity was first hinted at by Julia long ago 
[50, 51]). Note that the conjecture here is that the continuous group -Eio(lR) 
be a hidden symmetry group of classical supergravity. At the quantum level, 
i.e. for M-theory, one expects only a discrete version of E w , say E 10 (Z), 
to be a quantum symmetry. See [16] for recent work trying to extend the 
identification of [28] between roots of E w and symmetries of supergravity/M- 
theory beyond height 30, and for references about previous suggestions of a 
possible role of E w . For earlier appearances of the Weyl groups of the E 
series in the context of [/-duality see [61, 65, 3]. A series of recent papers 
[71, 68, 69, 35, 36] has also explored the possible role of E u (a non-hyperbolic 
extension of E w ) as hidden symmetry of M-theory. See also [41] for another 
approach to the search for symmetries of M theory. 

It is also tempting to assume that the Kac-Moody groups underlying 
the other (special) Einstein-matter systems discussed above might be hidden 
(solution-generating) symmetries. For instance, in the case of pure Einstein 
gravity in D = 4 spacetime, the conjecture is that AE 3 is such a symmetry 
of Einstein gravity. [This case, and the correspondence between the field 
variables and the coset ones is further discussed in [29].] 

To end this introductory summary, it is important to add a significant 
mathematical caveat. Most of the results mentioned above have been ob- 
tained by "physicists' methods", i.e. non rigorous arguments. [See, however, 
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the mathematical results of [67] concerning the dynamics of the Bianchi type 
IX homogeneous model.] The only exception concerns the "non chaotic", 
monotonic, Kasner-like behaviour of the "subcritical" systems. The first 
mathematical proof that a general solution of the four-dimensional Einstein- 
dilaton system has a non chaotic, Kasner-like behaviour was obtained by 
Andersson and Rendall [2]. By using, and slightly extending, the tools of 
this proof (Theorem 3 of [2], concerning certain Fuchsian systems), Ref. [30] 
has then given mathematical proofs of the Kasner-like behaviour of more gen- 
eral classes of Einstein-matter systems, notably the D-dimensional coupled 
Einstein-dilaton-p-form system when the dilaton couplings of the p-forms be- 
long to the subcritical domain mentioned above. Ref. [30] also gave a proof of 
the Kasner-like behaviour of pure gravity (Ricci tensor = 0) when the space- 
time dimension D > 11. For other mathematical studies on "non-chaotic" 
cases see Ref. [44]. 

An important mathematical challenge is to convert the physicists' ar- 
guments (summarized in the rest of this review) concerning the "chaotic 
billiard" structure of critical (and over critical) systems into precise mathe- 
matical statements. The recent work [28, 29], which provides a simplified 
derivation of the billiard picture and a rather clean decomposition of the field 
variables into "ODE-described chaotic" and "asymptotically frozen" pieces, 
might furnish a useful starting point for formulating precise mathematical 
statements. In this respect, let us note two things: (i) a good news is that 
numerical simulations have fully confirmed the BKL chaotic-billiard picture 
in several models [12, 11], (ii) a bad news is that the physicists' arguments 
might oversimplify the picture by neglecting "non generic" spatial points 
where some of the leading walls disappear, because the spatially-dependent 
coefficients measuring the strength of these walls happen to vanish. A simi- 
lar subtlety (vanishing at exceptional spatial points of wall coefficients) takes 
place in the subcritical, non chaotic case (where a finite number of collisions 
on the billiard walls is expected to channel the billiard ball in the "good", 
Kasner-like directions). See discussions of this phenomenon in [11, 66]. 

In the rest of this review, we shall outline the various arguments leading 
to the conclusions summarized above. For a more complete derivation of the 
billiard results see [29]. 
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2 Models and Gauge Conditions 



We consider models of the general form 



% M iv,<M {p) ] = J d D x 



R(g) - d. 



M<-. 



1 1 r K<t> F (P) p(p)Mi-M p +i 



V 



+ ... (2.1) 



where units are chosen such that \^G N = 1 (where G N is Newton's con- 
stant) and the spacetime dimension D = d + 1 is left unspecified. Besides 
the standard Einstein-Hilbert term the above Lagrangian contains a dilaton 
field and a number of p-form fields A M ^,,, Mp (for p > 0). For simplicity, 
we consider the case where there is only one dilaton, i.e. n = 1 in the no- 
tation of the Introduction. As a convenient common formulation we adopt 
the Einstein conformal frame and normalize the kinetic term of the dilaton 
(f> with weight one with respect to the Ricci scalar. The Einstein metric gMN 
has Lorentz signature ( — !-•••+) and is used to lower or raise the indices; 
its determinant is denoted by g. The p-form field strengths = dA^ are 
normalized as 

F mI-m p+1 = (j> + l ) d [MiA^ 2 ... Mv+i] = d Ml A { $ 2 ... Mp+i ±p permutations . (2.2) 

The dots in the action (2.1) indicate possible modifications of the field 
strength by additional Yang-Mills or Chapline-Manton-type couplings [13, 
17], such as F c = dC® - C^dB^ for two 2-forms and and a 
0-form as they occur in type IIB supergravity. Further modifications 

include Chern-Simons terms, as in the action for D = 11 supergravity [22]. 
The real parameter \ p measures the strength of the coupling of A^ to the 
dilaton. When p = 0, we assume that A ^ so that there is only one dilaton. 
This is done mostly for notational convenience. If there were other dilatons 
among the 0-forms, these should be separated off from the p-forms because 
they play a distinct role. They would define additional spacelike directions 
in the space of the (logarithmic) scale factors and would correspondingly 
increase the dimension of the relevant hyperbolic billiard. 

The metric gMN, the dilaton field(s) and the p-form fields A M \,,, Mp are 
a priori arbitrary functions of both space and time, on which no symme- 
try conditions are imposed. Nevertheless it will turn out that the evolution 
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equations near the singularity will be asymptotically the same as those of cer- 
tain homogeneous cosmological models. It is important to keep in mind that 
this simplification does not follow from imposing extra dimensional reduction 
conditions but emerges as a direct consequence of the general dynamics. 

Our analysis applies both to past and future singularities, and in particu- 
lar to Schwarzschild-type singularities inside black holes. To follow historical 
usage, we shall assume for definiteness that the spacelike singularity lies in 
the past, at finite distance in proper time. More specifically, we shall adopt 
a space-time slicing such that the singularity "occurs" on a constant time 
slice (t — in proper time). The slicing is built by use of pseudo- Gaussian 
coordinates defined by vanishing lapse N l = 0, with metric 

ds 2 = -{N{x°, x l )dx ) 2 + 9lj {x°, x l )dx l dxK (2.3) 

For simplicity, we shall work in the following with a coordinate spatial co- 
frame dx l , as indicated in Eq. (2.3). However, most of what we say would go 
through if we were working in a general time-independent spatial co-frame 
u l (x k ) = u>j(x k )dx 3 . [The usefulness of using a non-coordinate spatial co- 
frame will show up in the E w case discussed below.] In order to simplify 
various formulas later, we shall find it useful to introduce a rescaled lapse 
function 

N = N/^ (2.4) 

where g = det g^. [Note the distinction between spacetime quantities gMN, g 
and space ones g^, g. In the following, we work only with the spatial metric 
gij] We shall see that a useful gauge, within the Hamiltonian approach, is 
that defined by requiring 

N = p\ (2.5) 

where p 2 is a quadratic combination of the logarithms of the scale factors 
and the dilaton(s), which we will define below in terms of the Iwasawa de- 
composition of g^. 

Throughout this paper, we will reserve the label t for the proper time 

dt = -Ndx° = -N^gdx , (2.6) 

whereas the time coordinate associated with the special gauge Eq.(2.5) will 
be designated by T, viz. 

dT (2.7) 

p 2 s/g 



n 



Sometimes, it will also be useful to introduce the "intermediate" time coor- 
dinate r that would correspond to the gauge condition N — 1. It is explicitly 
defined by: 



At the singularity the proper time t is assumed to remain finite and to de- 
crease toward + . By contrast, the coordinates T and r both increase toward 
+00, as ensured by the minus sign in (2.6). Irrespective of the choice of co- 
ordinates, the spatial volume density g is assumed to collapse to zero at each 
spatial point in this limit. 

As for the p-form fields, we shall assume, throughout this paper, a gen- 
eralized temporal gauge, viz. 



where small Latin letter denote spatial indices from now on. This 

choice leaves the freedom of performing time-independent gauge transforma- 
tions, and therefore fixes the gauge only partially. 

3 Asymptotic dynamics in the general case 

The main ingredients of the derivation of the billiard picture are: 

• use of the Arnowitt-Deser-Misner Hamiltonian formalism [1] , 

• Iwasawa decomposition of the metric, i.e. — > (f3 a ,Af a i), 

• decomposition of f3^ = (f3 a , 0) into radial (p) and angular (7^) parts, 



• use of the pseudo-Gaussian gauge (2.5), i.e. of the time coordinate T 
as the evolution parameter. 

More explicitly, with the conventions already described before, we as- 
sume that in some spacetime patch, the metric is given by (2.3) (pseudo- 
Gaussian gauge), such that the local volume g collapses at each spatial point 
as x° — > +00, in such a way that the proper time t tends to + . We work in 
the Arnowitt-Deser-Misner Hamiltonian formalism [1], i.e. with first order 



dr = = p 2 dT. 

^9 



(2.8) 




(2.9) 



and 
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evolution equations in the phase-space of the system. For instance, the grav- 
itational degrees of freedom are initially described by the metric and its 
conjugate momentum 7r u . We systematically use the Iwasawa decomposition 
of the spatial metric g^. Namely, we write 

g %3 = Y^e-^ a M a t M a 3 (3.1) 

a=l 

where Af is an upper triangular matrix with l's on the diagonal, and where 
we recall that d = D — 1 denotes the spatial dimension. We shall also need 
the Iwasawa coframe {6 a } 

Q a =M a i dx i (3.2) 
as well as the vectorial frame {e a } dual to the coframe {9 a }: 

e a =A/\A (3.3) 

where the matrix M l a is the inverse of N a i, i.e., A/" a jA/\ = It is again an 
upper triangular matrix with l's on the diagonal. 

The Iwasawa decomposition allows us to replace the d(d+ l)/2 variables 
gtj by the d+d(d-l)/2 variables (f3 a , M a i). Note that (f3 a , N a i) are ultralocal 
functions of g^, that is they depend, at each spacetime point, only on the 
value of gij at that point, not on its derivatives. This would not have been 
the case if we had used a "Kasner frame" (i.e. a frame diagonalizing -k^ with 
respect to g^) instead of an Iwasawa one. The transformation g — > ()3,Af) 
then defines a corresponding transformation of the conjugate momenta, as 
we will explain below. We then augment the definition of the /3's by adding 
the dilaton field, i.e. (3^ = (/3 a , 0). The d + n dimensional space of the /5 M 's 
(where, as we said, we consider only one dilaton, i.e. n = 1) comes equipped 
with a canonical Lorentzian metric G^ u defined by 

d d 

da 2 = G^ v df3^df3 u = Y^W 1 ? ~ E + ( 3 - 4 ) 

a=l a=l 

This (flat) Lorentzian metric in the auxiliary /9-space (of signature — K . • +) 
plays an essential role in our study. Note that in models without dilaton 

3 If we were working with a non-coordinate spatial co-frame u l = uj^(x k )dx^ , we would 
use the Iwasawa decomposition of the frame components of the spatial metric: di 2 — 
g lJ uJ l uJ 1 . 
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(such as pure gravity in spacetime dimension d+ 1) one has a rf-dimensional 
Lorentzian space, with (/S' 1 ) = ((3 a ) and metric G iW dfi ,x dfi v = ^2 a {d(3 a ) 2 — 

We then decompose f3^ in hyperbolic polar coordinates (p, 7^), z.e 

0" = P7 / \ (3.5) 

where 7^ are coordinates on the future sheet of the unit hyperboloid 4 , which 
are constrained by 

G^7 M 7^7^ = -1 (3.6) 
and p is the timelike variable 

p 2 = -G, u p^P" = -p^,>0. (3.7) 

This decomposition naturally introduces the unit hyperboloid ("7-space"), 
see Fig. 1, which is a realization of the m-dimensional hyperbolic (Lobachev- 
skii) space H m , with m = d — 1 + n if there are n > dilatons. 

In terms of the "polar" coordinates p and 7^, the metric in /5-space be- 
comes 

da 2 = -dp 2 + p 2 dT 2 (3.8) 

where dTi 2 is the metric on the 7-space H m . 

Note that p is also an ultralocal function of the configuration variables 
(gij,<f>)- We assume that the hyperbolic coordinate p can be used every- 
where in a given region of space near the singularity as a well-defined (real) 
quantity which tends to +00 as we approach the singularity. We then de- 
fine the slicing of spacetime by imposing the gauge condition (2.5). The 
time coordinate corresponding to this gauge is called T as above (see (2.7)). 
Our aim is to study the asymptotic behaviour of all the dynamical variables 

(3{T),N{T), as T — > +00 (recall that this limit also corresponds to t — > 0, 

^J~g — * 0, p — > +00, with (3^ going to infinity inside the future light cone). 
Of course, we must also ascertain the self-consistency of this limit, which we 
shall refer to as the "BKL limit" . 

3.1 Arnowitt-Deser-Misner Hamiltonian action 

To focus on the features relevant to the billiard picture, we assume here that 
there are no Chern-Simons terms or couplings of the exterior form gauge 

4 Indeed, one finds that, near a spacelike singularity /3 M tends to future timelike infinity. 
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fields through a modification of the curvatures F&\ which are thus taken to 
be Abelian, = dA^\ See [29] for a proof that these interaction terms 
do not change the analysis. The Arnowitt-Deser-Misner Hamiltonian action 
in any pseudo-Gaussian gauge, and in the temporal gauge (2.9), reads [1] 

J dx° J d d x ^% • ^ • i £ vrg)"^.,, - Hj (3.9) 

where the Hamiltonian density H is 

H = NH (3.10) 
H = fC + M (3.11) 

( 3 - 12 ) 



2p! 

p 



where i? is the spatial curvature scalar. The dynamical equations of motion 
are obtained by varying the above action with respect to the spatial metric 
components, the dilaton, the spatial p-form components and their conjugate 
momenta. In addition, there are constraints on the dynamical variables, 

H « ("Hamiltonian constraint"), (3.14) 
Hi ~ ("momentum constraint"), (3.15) 
</^)" Jp_1 ~ ("Gauss law" for each p-form) (3.16) 



with 



T~ti 



v 



where the subscript \j stands for the spatially covariant derivative. 
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Let us now see how the Hamiltonian action gets transformed when one 
performs, at each spatial point, the Iwasawa decomposition (3.1) of the spa- 
tial metric. The kinetic terms of the metric and of the dilaton in the La- 
grangian (2.1) are given by the quadratic form 

X>/? a ) 2 - fl>/3 a ) + d(j> 2 + l - Y^eW-^ (dAf^M^) 2 (3.19) 

a=l \a=l / a<b 

where we recall that N l a denotes, as in (3.3), the inverse of the triangular ma- 
trix M a i appearing in the Iwasawa decomposition (3.1) of the spatial metric 
gij. This change of variables corresponds to a point canonical transformation, 
which can be extended to the momenta in the standard way via 

7T% EE J2 + J2 Pt ^" ( 3 - 20 ) 

a a 

Note that the momenta 

OJ V i ^ 

conjugate to the non-constant off-diagonal Iwasawa components Af a i, are only 
defined for a < i; hence the second sum in (3.20) receives only contributions 
from a < i. 

We next split the Hamiltonian density 5 H (3.10) in two parts, one denoted 
by 7i , which is the kinetic term for the local scale factors f3^ (including 
dilatons) already encountered in section 3, and a "potential density" (of 
weight 2) denoted by V, which contains everything else. Our analysis below 
will show why it makes sense to group the kinetic terms of both the off- 
diagonal metric components and the p-forms with the usual potential terms, 
i.e. the term Ai in (3.10). [Remembering that, in a Kaluza-Klein reduction, 
the off-diagonal components of the metric in one dimension higher become 
a one-form, it is not surprising that it might be useful to group together the 
off-diagonal components and the p-forms.] Thus, we write 

H = H + V (3.22) 

5 We use the term "Hamiltonian density" to denote both H and H. Note that H is a 
usual spatial density (of weight 1, i.e. the same weight as y/g), while H = ^fgH/N is a 
density of weight 2 (like g = (y/g) 2 )- Note also that 7r y is of weight 1, while N = N/^/g 
is of weight — 1 . 
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with the kinetic term of the (3 variables 

Ho = ^%7r„ (3.23) 

where G^ u denotes the inverse of the metric of Eq. (3.4). In other words 
the right hand side of Eq. (3.23) is defined by 

where = (n a , n^) are the momenta conjugate to f3 a and 0, respectively, 
i.e. 

TV = 2N' 1 G lJl J l/ = *gJ^. (3.25) 
The total (weight 2) potential density, 

V = V 5 + V G + ]TV P + V , (3.26) 

v 

is naturally split into a "centrifugal" part linked to the kinetic energy of 
the off-diagonal components (the index "S"' referring to "symmetry", as dis- 
cussed below) 

V5 = ^e- 2 ^°)(P^) 2 , (3.27) 

a<b 

a "gravitational" (or "curvature") potential 

V G = -gR, (3.28) 

and a term from the p-forms, 

V(p) = Vg } + VST" (3.29) 
which is a sum of an "electric" and a "magnetic" contribution 

K) = -2^T 7r 5)" J, ' 7r (p)ii-ip (3-30) 

*>) _ 2(p + l)! 5 • 1 J 

Finally, there is a contribution to the potential linked to the spatial gradients 
of the dilaton: 

V = gg ij d4d j( t>. (3.32) 
We will analyze below in detail these contributions to the potential. 
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3.2 Appearance of sharp walls in the BKL limit 

In the decomposition of the Hamiltonian given above, we have split off the 
kinetic terms of the scale factors (3 a and of the dilaton (3 d+1 = from the 
other variables, and assigned the off-diagonal metric components and the p- 
form fields to various potentials, each of which is a complicated function of 
(3^,Af a i, P' l a , Aj^,,,j p , 7r^ " Jp and of some of their spatial gradients. The reason 
why this separation is useful is that, as we are going to show, in the BKL 
limit, and in the special Iwasawa decomposition which we have adopted, the 
asymptotic dynamics is governed by the scale factors whereas all other 
variables "freeze". Thus, in the asymptotic limit, we have schematically 

v(p^N\, P\, Afl Jp , Trg"^, . . . ) — VUY) (3.33) 

where V 00 {l fl ) stands for a sum of certain "sharp wall potentials" which 
depend only on the angular hyperbolic coordinates 7 M = (3^ / p. As a con- 
sequence, the asymptotic dynamics can be described as a "billiard" in the 
hyperbolic space of the 7 M 's, whose walls (or "cushions") are determined by 
the energy of the fields that are asymptotically frozen. 

This reduction of the complicated potential to a much simpler "effective 
potential" Voo(7 M ) follows essentially from the exponential dependence of 
V on the diagonal Iwasawa variables from its independence from the 
conjugate momenta of the /3's, and from the fact that the radial magnitude 
p of the /9's becomes infinitely large in the BKL limit. 

To see the essence of this reduction, with a minimum of technical com- 
plications, let us consider a potential density (of weight 2) of the general 
form 

V(/3, Q, P) = J2 c a(Q, P) exp ( - 2w A (P)) (3.34) 

A 

where (Q, P) denote the remaining phase space variables (that is, other than 
((3,77/3)). Here Wa{/3) = Wa^^ are certain linear forms which depend only 
on the (extended) scale factors, and whose precise form will be derived in 
the following section. Similarly we shall discuss below the explicit form of 
the pre-factors ca, which will be some complicated polynomial functions of 
the remaining fields, i.e. the off-diagonal components of the metric, the p- 
form fields and their respective conjugate momenta, and of some of their 
spatial gradients. The fact that the wa{(3) depend linearly on the scale 
factors is an important property of the models under consideration. A 
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second non-trivial fact is that, for the leading contributions, the pre-factors 
are always non-negative, i.e. c ^ ading > o. Since the values of the fields for 
which ca = constitute a set of measure zero, we will usually make the 
"genericity assumption" ca > for the leading terms in the potential V 6 . 
The third fact following from the detailed analysis of the walls that we shall 
exploit is that all the leading walls are timelike, i. e. their normal vectors (in 
the Minkowski /3-space) are spacelike. 

When parametrizing (3^ in terms of p and 7 M , or equivalently 

X = \np=\\n{-G^n (3.35) 

(with conjugate momentum tt\), and 7^, the part of the Hamiltonian des- 
cribing the kinetic energy of the /3's, H = NHq, takes the form 

H o = 7^ W + <] ■ ( 3 - 36 ) 

Choosing the gauge (2.5) to simplify the kinetic terms H we end up with an 
Hamiltonian of the form 

tf(A,7T A , 7 ,7r 7 ,Q,P) = NH (3.37) 
= \ [-Ta + <] + P 2 Yl c ^> P ) ex P ( - ^(7)) 

A 

where is the kinetic energy of a particle moving on H m . In (3.37) and 
below we shall regard A as a primary dynamical variable (so that p = e A ). 

The essential point is now that, in the BKL limit, A — > +00 i.e. p — > +00, 
each term p 2 exp ( — 2pw^(7)) becomes a sharp wall potential, i.e. a function 
of Wa(i) which is zero when ^(7) > 0, and +00 when 11^(7) < 0. To 
formalize this behaviour we define the sharp wall 9-function 7 as 

e(s):=|° l X< l (3.38) 
v ; \ +00 if x > v ; 

A basic formal property of this 6-function is its invariance under multiplica- 
tion by a positive quantity. With the above assumption checked below that 

6 As mentioned in the Introduction understanding the effects of the possible failure of 
this assumption is one of the subtle issues in establishing a rigorous proof of the BKL 
picture. 

7 Onc should more properly write 600(2:), but since this is the only step function en- 
countered in this article, we use the simpler notation Q(x). 
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all the relevant prefactors ca(Q,P) are positive near each leading wall, we 
can formally write 

lim \c A (Q,P)p 2 exp(-pw A (-f)} = c A (Q, P)G( - 2w A ( 1 )) = 6( - 2w A (>y)). 

(3.39) 

Of course, Q{—2w A { , y)) = 0(— 11^(7)), but we shall keep the extra factor of 2 
to recall that the arguments of the exponentials, from which the O-functions 
originate, come with a well-defined normalization. Therefore, the limiting 
Hamiltonian density reads 

#00 (A, vr A , 7, vr 7 , Q,P) = \ + tt*] + ^ ©( ~ 2u^(t)), (3.40) 

where A' runs over the dominant walls. The set of dominant walls is defined 
as the minimal set of wall forms which suffice to define the billiard table, i.e. 
such that the restricted set of inequalities {w A >(^) > 0} imply the full set 
{w A (^) > 0}. The concept of dominant wall will be illustrated below. [Note 
that the concept of "dominant" wall is a refinement of the distinction, which 
will also enter our discussion, between a leading wall and a subleading one.] 
The crucial point is that the limiting Hamiltonian (3.40) no longer de- 
pends on A, Q and P. Therefore the Hamiltonian equations of motion for 
A, Q and P tell us that the corresponding conjugate momenta, i.e. 7r A , P and 
Q, respectively, all become constants of the motion in the limit A — > +00. 
The total Hamiltonian density is also a constant of the motion (which 
must be set to zero). The variable A evolves according to dX/dT = — jir\. 
Hence, in the limit, A is a linear function of T. The only non-trivial dynamics 
resides in the evolution of (7, 7r 7 ) which reduces to the sum of a free (non 
relativistic) kinetic term vr^/4 and a sum of sharp wall potentials, such that 
the resulting motion of the 7's indeed constitutes a billiard, with geodesic 
motion on the unit hyperboloid H m interrupted by reflections on the walls 
defined by w A {j) = 0. These walls are hyperplanes (in the sense of hyper- 
bolic geometry) because they are geometrically given by the intersection of 
the unit hyperboloid (3^(3^ = — 1 with the usual Minkowskian hyperplanes 
w A (P) = 0. 

See Appendix A of [29] for the study of a toy model which shows in 
more detail how the asymptotic constancy of the "off diagonal" phase-space 
variables (Q, P) arises. 

Geodesic motion in a billiard in hyperbolic space has been much studied. 
It is known that this motion is chaotic or non-chaotic depending on whether 
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the billiard has finite or infinite volume [62, 42, 72, 37]. In the finite volume 
case, the generic evolution exhibits an infinite number of collisions with the 
walls with strong chaotic features ("oscillating behavior"). 

By contrast, if the billiard has infinite volume, the evolution is non- 
chaotic. For a generic evolution, there are only finitely many collisions with 
the walls. The system generically settles after a finite time in a Kasner-like 
motion that lasts all the way to the singularity, see Fig. 1. 

The above derivation relied on the use of hyperbolic polar coordinates 
(p, 7). This use is technically useful in that it represents the walls as be- 
ing located at an asymptotically fixed position in hyperbolic space, namely 
wa{i) = 0. However, once one has derived the final result (3.40), one 
can reexpress it in terms of the original variables (3^, which run over a lin- 
ear (Minkowski) space. Owing to the linearity of wa{/3) = pwa(i) in this 
Minkowskian picture, the asymptotic motion takes place in a "polywedge", 
bounded by the hyperplanes wa(/3) = 0. The billiard motion then consists 
of free motions of /3 M on straight lightlike lines within this polywedge, which 
are interrupted by specular reflections off the walls. [See formula (5.3) below 
for the explicit effect of these reflections on the components of the velocity 
vector of the /5-particle.] Indeed, when going back to /5-space (i.e. before 
taking the BKL limit), the dynamics of the scale factors at each point of 
space is given by the Hamiltonian 



H({3», tt m ) =NH = N 



\Cr^ v + ^ c A exp ( - 2w A ((3)) 



(3.41) 



The /9-space dynamics simplifies in the gauge N — 1, corresponding to the 
time coordinate r. In the BKL limit, the Hamiltonian (3.41) takes the lim- 
iting form (in the gauge N — 1) 

HocW, tt m ) = \G^n u + ^ 8( - 2 «M/3)) (3-42) 

A' 

where the sum is again only over the dominant walls. When taking "equal 
time" slices of this polywedge (e.g., slices on which Ej/3 l is constant), it is 
clear that with increasing time (i.e increasing Ej/?, or increasing r or p) the 
walls recede from the observer. The /3-space picture is useful for simplifying 
the mathematical representation of the dynamics of the scale factors which 
takes place in a linear space. However, it is inconvenient both for proving 
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that the exponential walls of (3.41) do reduce, in the large limit to 
sharp walls, and for dealing with the dynamics of the other phase-space 
variables (Q, P), whose appearance in the coefficients ca has been suppressed 
in (3.41) above. Let us only mention that, in order to prove, in this picture, 
the freezing of the phase space variables (Q, P) one must consider in detail 
the accumulation of the "redshifts" of the energy-momentum 7r M of the (3- 
particle when it undergoes reflections on the receding walls, and the effect 
of the resulting decrease of the magnitudes of the components of 7r M on the 
evolution of (Q,P). By contrast, the 7-space picture that we used above 
allows for a more streamlined treatment of the effect of the limit p — > +00 
on the sharpening of the walls and on the freezing of (Q,P). 

In summary, the dynamics simplifies enormously in the asymptotic limit. 
It becomes ultralocal in that it reduces to a continuous superposition of 
evolution systems (depending only on a time parameter) for the scale factors 
and the dilatons, at each spatial point, with asymptotic freezing of the off- 
diagonal and p-form variables. This ultralocal description of the dynamics 
is valid only asymptotically. It would make no sense to speak of a billiard 
motion prior to this limit, because one cannot replace the exponentials by 
O-functions. Prior to this limit, the evolution system for the scale factors 
involves the coefficients in front of the exponential terms, and the evolution 
of these coefficients depends on various spatial gradients of the other degrees 
of freedom. However, one may contemplate setting up an expansion in which 
the sharp wall model is replaced by a model with exponential ( "Toda-like" ) 
potentials, and where the evolution of the quantities entering the coefficients 
of these "Toda walls" is treated as a next to leading effect. See [28] and below 
for the definition of the first steps of such an expansion scheme for maximal 
supergravity in eleven dimensions. 

3.3 Constraints 

We have just seen that in the BKL limit, the evolution equations become 
ordinary differential equations with respect to time. Although the spatial 
points are decoupled in the evolution equations, they are, however, still cou- 
pled via the constraints. These constraints just restrict the initial data and 
need only be imposed at one time, since they are preserved by the dynamical 
equations of motion. Indeed, one easily finds that, in the BKL limit, 

H = (3.43) 
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since [H(x),7i(x')] = in the ultralocal limit. This corresponds simply to 
the fact that the collisions preserve the lightlike character of the velocity 
vector. Furthermore, the gauge constraints (3.16) are also preserved in time 
since the Hamiltonian constraint is gauge- invariant. In the BKL limit, the 
momentum constraint fulfills 

n k {x) = d k n « o. (3.44) 

It is important that the restrictions on the initial data do not bring dangerous 
constraints on the coefficients of the walls in the sense that these may all take 
non-zero values. For instance, it is well known that it is consistent with the 
Gauss law to take non-vanishing electric and magnetic energy densities; thus 
the coefficients of the electric and magnetic walls are indeed generically non- 
vanishing. In fact, the constraints are essentially conditions on the spatial 
gradients of the variables entering the wall coefficients, not on these variables 
themselves. In some non-generic contexts, however, the constraints could 
force some of the wall coefficients to be zero; the corresponding walls would 
thus be absent. [E.g., for vacuum gravity in four dimensions, the momentum 
constraints for some Bianchi homogeneous models force some symmetry wall 
coefficients to vanish. But this is peculiar to the homogeneous case.] 

It is easy to see that the number of arbitrary physical functions involved 
in the solution of the asymptotic BKL equations of motion is the same as 
in the general solution of the complete Einstein-matter equations. Indeed, 
the number of constraints on the initial data and the residual gauge freedom 
are the same in both cases. Further discussion of the constraints in the BKL 
context may be found in [2, 30]. 

3.4 Consistency of BKL behaviour in spite of the in- 
crease of spatial gradients 

The essential assumption in the BKL analysis is the asymptotic dominance of 
time derivatives with respect to space derivatives near a spacelike singularity. 
This assumption has been mathematically justified, in a rigorous manner, in 
the cases where the billiard is of infinite volume, i.e. in the (simple) cases 
where the asymptotic behaviour is not chaotic, but is monotonically Kasner- 
like [2, 30]. 

On the other hand, one might a priori worry that this assumption is 
self-contradictory in those cases where the billiard is of finite volume, when 
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the asymptotic behaviour is chaotic, with an infinite number of oscillations. 
Indeed, it has been pointed out [55, 6] that the independence of the billiard 
evolution at each spatial point will have the effect of infinitely increasing 
the spatial gradients of various quantities, notably of the local values of the 
Kasner exponents p^(x). This increase of spatial gradients towards the singu- 
larity has been described as a kind of turbulent behaviour of the gravitational 
field, in which energy is pumped into shorter and shorter length scales [55, 6], 
and, if it were too violent, it would certainly work against the validity of the 
BKL assumption of asymptotic dominance of time derivatives. For instance, 
in our analysis of gravitational walls in the following section, we will en- 
counter subleading walls, whose prefactors depend on spatial gradients of 
the logarithmic scale factors (3. 

To address the question of consistency of the BKL assumption we need to 
know how fast the spatial gradients of /3, and of similar quantities entering 
the prefactors, grow near the singularity. We refer to [29] for a discussion 
of this issue. The result is that the chaotic character of the billiard indeed 
implies an unlimited growth of the spatial gradients of (3, but that this growth 
is only of polynomial order in p 

8,(3 = 0(l)p 2 , (3.45) 

where the coefficient 0(1) is a chaotically oscillating quantity. This poly- 
nomial growth of di/3 (and of its second-order spatial derivatives) entails a 
polynomial growth of the prefactors of the sub-dominant walls. Because it 
is polynomial (in p), this growth is, however, negligible compared to the 
exponential (in p) behaviour of the various potential terms. It does not 
jeopardize our reasoning based on keeping track of the various exponential 
behaviours. As we will see the potentially dangerously growing terms that 
we have controlled here appear only in subdominant walls. The reasoning 
of the Appendix of [29] shows that the prefactors of the dominant walls are 
self-consistently predicted to evolve very little near the singularity. 

We conclude that the unlimited growth of some of the spatial gradients 
does not affect the consistency of the BKL analysis done here. This does 
not mean, however, that it will be easy to promote our analysis to a rigorous 
mathematical proof. The main obstacle to such a proof appears to be the 
existence of exceptional points, where a prefactor of a dominant wall happens 
to vanish, or points where a subdominant wall happens to be comparable to 
a dominant one. Though the set of such exceptional points is (generically) 
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of measure zero, their density might increase near the singularity because of 
the increasing spatial gradients. This situation might be compared to the 
KAM (Kolmogorov-Arnold-Moser) one, where the "bad" tori have a small 
measure, but are interspersed densely among the "good" ones. 

4 Walls 

The decomposition (3.22) of the Hamiltonian gives rise to different types 
of walls, which we now discuss in turn. Specifically, we will derive explicit 
formulas for the linear forms wa(/3) and the field dependence of the pre- 
factors ca entering the various potentials. 

4.1 Symmetry walls 

We start by analyzing the effects of the off-diagonal metric components which 
will give rise to the so-called "symmetry walls" . As they originate from the 
gravitational action they are always present. The relevant contributions to 
the potential is the centrifugal potential (3.27). When comparing (3.27) to 
the general form (3.34) analyzed above, we see that firstly the summation 
index A must be interpreted as a double index (a, b), with the restriction 
a < b, secondly that the corresponding prefactor is c a b = (P\-^ b j) 2 is auto- 
matically non-negative (in accordance with our genericity assumptions, we 
shall assume c ab > 0). The centrifugal wall forms read: 

= = ~ P a (a < b). (4.1) 

We refer to these wall forms as the "symmetry walls" for the following reason. 
When applying the general collision law (5.3) derived below to the case of the 
collision on the wall (4.1) one easily finds that its effect on the components of 
the velocity vector v M is simply to permute the components v a and v b , while 
leaving unchanged the other components fi ^ a,b. 

The hyperplanes w^ ab ^(/3) = (i.e. the symmetry walls) are timelike since 

G^wf ab)ll w s {ab)v = +2. (4.2) 

This ensures that the symmetry walls intersect the hyperboloid G^ v l3^f3 v = 
— 1, 'Yl ia f3 a > 0. The symmetry billiard (in /9-space) is defined to be the 
region of Minkowski space determined by the inequalities 

>0, (4.3) 
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with J2aft a > (*- e - by the region of /3-space where the 9 functions are 
zero). Its projection on the hyperbolic space H m is defined by the inequalities 



The explicit expressions above of the symmetry wall forms also allow us 
to illustrate the notion of a "dominant wall" defined above. Indeed, the 
d(d — l)/2 inequalities (4.3) already follow from the following minimal set of 
d — 1 inequalities 



More precisely, each linear form which must be positive in (4.3) can be written 
as a linear combination, with positive (in fact, integer) coefficients of the 
linear forms entering the subset (4.4). For instance, P 3 ~P 1 = (/3 3 —j3 2 ) + (f3 2 — 
P 1 ), etc. As discussed in [29] this result can be reinterpreted by identifying 
the dominant linear forms entering (4.4) with the simple roots of SL{n, M). 

4.2 Curvature (gravitational) walls 

Next we analyze the gravitational potential, which requires a computation 
of curvature. To that end, one must explicitly express the spatial curvature 
in terms of the scale factors and the off-diagonal variables M a i. Again, the 
calculation is most easily done in the Iwasawa frame (3.2). We use the short- 
hand notation A a = e _/3 " for the (Iwasawa) scale factors. Let C a hc (x) be the 
structure functions of the Iwasawa basis {0 a }, viz. 



where d is the spatial exterior differential. The structure functions obviously 
depend only on the off-diagonal components M a i) but not on the scale factors. 
Using the Cartan formulas for the connection one-form u a b , 



«f«b)(7) > 0. 



P 2 - P 1 > 0, P 3 - P 2 > 0, • • • , (3 d - p d ~ l > 0. 



(4.4) 



= - l -C\ c Q h NQ c 



(4.5) 




(4.6) 



b 



d^ab = U ab + U bt 



HI 



(4.7) 



where uj ab = ^ a c^ c b, and 



lab 



= S ab A 2 a = 



exp(-2p a )S ab 



(4.8) 
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is the metric in the frame {9 a }, one finds 



, , c 1 / n b A-l , A d „ c \ , 

+ E 2^2 [M^c),* + M^c),<* - M^),c] (4.9) 

In the last bracket above, the commas denote the frame derivatives d a = 
A/\<9,. 

The Riemann tensor R c de f, the Ricci tensor R de and the scalar curvature 
R are obtained through 

n\ = ^ + 5> a c Ao/ fe (4.10) 

c 

= \Y.R\e f o e ^ef (4.ii) 

e,f 

where Q a b is the curvature 2-form and 

Rab = ^2 RC acbi R = ^ ~T2 Raa - ( 4 - 12 ) 

c a a 

Direct, but somewhat cumbersome, computations yield 

* = -z E -A^) 2 + E -^(^ ^. 5C > <?) ( 4 - 13 ) 

a,b,c b c a a 

where F a is some complicated function of its arguments whose explicit form 
will not be needed here. The only property of F a that will be of importance is 
that it is a polynomial of degree two in the derivatives df3 and of degree one 
in d 2 f3. Thus, the exponential dependence on the /3's which determines the 
asymptotic behaviour in the BKL limit, occurs only through the A^-terms 
written explicitly in (4.13). 

In (4.13) one obviously has b ^ c because the structure functions C a bc are 
antisymmetric in the pair [be] . In addition to this restriction, we can assume, 
without loss of generality, that a ^ b, c in the first sum on the right-hand side 
of (4.13). Indeed, the terms with either a = b or a = c can be absorbed into 
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a redefinition of F a . We can thus write the gravitational potential density 
(of weight 2) as 

V G = -gR = -J^e- 2a ^\C\ c f - e~ 2 ^ {f3) F a (4.14) 

a,b,c a 

where the prime on indicates that the sum is to be performed only over 
unequal indices, i.e. a ^ b,b ^ c,c ^ a, and where the linear forms a abc {(3) 
and fJ> a (/3) are given by 

a abc {{3) = 2f3 a + F (a^b, b^c, c^a) (4.15) 

e^a,b,c 

and 

respectively. Note that a abc is symmetric under the exchange of b with c, but 
that the index a plays a special role. 

Comparing the result (4.14) to the general form (3.34) we see that there 
are, a priori, two types of gravitational walls: the a-type and the /x-type. It 
is shown in [29] that the a-walls dominate the /i-walls and, more precisely 
that (J, a (/3) > within the entire future light cone of the /3's. [The proof uses 
the fact that each linear form /x a (/3) is lightlike.] 

From these considerations we deduce the additional constraints 

aabc(P) > (D > 3) (4.17) 

besides the symmetry inequalities (4.4). The hyperplanes a abc ((3) = are 
called the "curvature" or "gravitational" walls. Like the symmetry walls, 
they are timelike since 

G^(a abc ) ll (a abc ) u = +2. (4.18) 

The restriction D > 3 is due to the fact that in D = 3 spacetime dimensions, 
the gravitational walls a abc {l3) = are absent, simply because one cannot find 
three distinct spatial indices. This is, of course, in agreement with expec- 
tations, because gravity in three spacetime dimensions has no propagating 
degrees of freedom (gravitational waves). 
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4.3 p-form walls 

While none of the wall forms considered so far involved the dilatons, the 
electric and magnetic ones do as we shall now show. To make the notation 
less cumbersome we will omit the super-(or sub-)script (p) on the p-form 
fields in this subsection. 

4.3.1 Electric walls 

The electric potential density can be written as 

V t>) = —\ E e- 2e ^- a ^\£ ai - ap ) 2 (4.19) 

^ ai ,0,2 ,■ ■ ■ ,a p 

where S ai "' ap are the components of the electric field ft 11 '" 1 ? in the basis {6 a } 

£ ai - a " = M ai h N a2 h ■ ■ ■H ap j y i '" jp (4.20) 

(recall our summation conventions for spatial cordinate indices) and where 
e a i-a p (P) are the electric wall forms 8 

e ai ... ap ((3) = (3 a > + ■ ■ ■ + (3 a * + ^cj>. (4.21) 

Here the indices a/s are all distinct because S ai "' ap is completely antisym- 
metric. The variables S ai "' ap do not depend on the f3^. It is thus rather easy 
to take the BKL limit. The exponentials in (4.19) are multiplied by positive 
factors which generically are different from zero. Thus, in the BKL limit, 
V'fp) becomes 

K)^ E ©[-2e ai ... ap (/3)]. (4.22) 

ai<a,2<-<a p 

The transformation from the variables (A/" a j, P l a , Aj 1 ...j p , 7r ?1 '" J,p ) to the vari- 
ables (Af a i, V\, A ai -a v -> £ ai "' ap ) is a point canonical transformation whose 
explicit form is obtained from 

+ E y 1 ^^ = E^ a * + E E V ai - ap A,.. ap . 

a P a P ai,...,a p 

(4.23) 

8 A sign error in the last term cx A p crept in the corresponding equation of [29]. 
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The new momenta V' l a conjugate to M a i differ from the old ones P l a by 
terms involving £, M and A since the components A ai ... ap of the p-forms in 
the basis {9 a } depend on the A/"'s, 

A = Af jl ■ ■ ■ Af jp A ■ ■ 

However, it is easy to see that these extra terms do not affect the symmetry 
walls in the BKL limit. 



4.3.2 Magnetic walls 

The magnetic potential is dealt with similarly. Expressing it in the {0 a }- 
frame, one obtains 

VZT= „ * J2 ^ 2mai ^ +l(P \Fa 1 ...a p+1 ? (4.24) 



(P) 



2(p+l)! 

ai,02,-")O p +i 



where J 7 ai -a p+1 are the components of the magnetic field F mi ... m +1 in the 
basis {9 a }, 

F ai -a p+1 = -AP'V ■ ■■N^\ p+1 F n ... 3p+1 . (4.25) 
The m ai ... ap+1 (P) are the magnetic linear forms 9 



L ai---a p j r i 

b£{ai,a2,---a p +i} 



where again all a/s are distinct. One sometimes rewrites m ai ... ap+1 ((3) as 
m ap+2 ... ad , where {a p+2 ,a p+3 , ■ ■ ■ ,a d } is the set complementary to 
{a 1: a 2 , ■ ■ -a p+1 }; e.g., 

A 



m 12 .-d- 



P -i = P 1 + --- + P^ 1 - -i = m d _ p ... d . (4.27) 



Of course, the components of the exterior derivative T of ^4 in the non- 
holonomic frame {8 a } involves the structure coefficients, i.e. J 7 ai ---a p+1 = 
d[ ai A a2 ... ap+1 ] + C„4-terms where d a = M l a di is the frame derivative. 

Again, the BKL limit is quite simple and yields (assuming generic mag- 
netic fields) 

K) 9 ^ E e[-2m ai ... ad _ p _ 1 (/3)]. (4.28) 

ai<---<ad- p -i 

9 A sign error in the last term oc A p crept in the corresponding equation of [29]. 
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Just as the off-diagonal variables, the electric and magnetic fields freeze in the 
BKL limit since the Hamiltonian no longer depends on the p-form variables. 
These drop out because one can rescale the coefficient of any G-function to 
be one (when it is not zero), thereby absorbing the dependence on the p-form 
variables. 

The scale factors are therefore constrained by the further "billiard" con- 
ditions 

e ai -..aM > 0, ^...^(,3) >0. (4.29) 

The hyperplanes e ai ... ap (/3) = and m ai ... ad _ p _ 1 (j3) = are called "electric" 
and "magnetic" walls, respectively. Both walls are timelike because their 
gradients are spacelike, with squared norm 



p(d — p — 1 ) /A 



d- 1 



2 



+ (?) >0. (4.30) 



(For D — 11 supergravity, we have d — 10, p — 3 and X p = and thus the 
norm is equal to +2.) This equality explicitly shows the invariance of the 
norms of the p-form walls under electric-magnetic duality. 



4.4 Subleading walls 

The fact that the leading walls originating from the centrifugal, gravitational 
and p-form or are all timelike, is an important ingredient of the overall BKL 
picture. We refer to [29] for a discussion of the other contributions to the 
Hamiltonian (including Chern-Simons and Yang-Mills couplings) and for a 
proof that they contribute only subleading walls. 



5 Einstein (or cosmological) billiards 

Let us summarize the findings above. The dynamics in the vicinity of a 
spacelike singularity is governed by the scale factors, while the other vari- 
ables (off-diagonal metric components, p-form fields) tend to become mere 
"spectators" which get asymptotically frozen. This simple result is most eas- 
ily derived in terms of the hyperbolic polar coordinates (p, 7), and in the 
gauge (2.5). In this picture, the essential dynamics is carried by the angular 
variables 7 which move on a fixed billiard table, with cushions defined by the 
dominant walls wa>{^)- One can refer to this billiard as being an "Einstein 
billiard" (or a "cosmological billiard"). 
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It is often geometrically more illuminating to "unproject" the billiard 
motion in the full Minkowski space of the extended scale factors (3^. In that 
picture, the asymptotic evolution of the scale factors at each spatial point 
reduces to a zigzag of null straight lines with respect to the metric G^dft^dfl" . 
The straight segments of this motion are interrupted by collisions against the 
sharp walls 

w A {(3) = w A ^ = (5.1) 

defined by the symmetry, gravitational and p-form potentials, respectively. 
As we showed all these walls are timelike, i.e. they have spacelike gradients: 

G^WApWAv > 0. (5.2) 

Indeed, the gradients of the symmetry and gravitational wall forms have 
squared norm equal to +2, independently of the dimension d. By contrast, 
the norms of the electric and magnetic gradients, which are likewise positive, 
depend on the model. As we saw, there also exist subdominant walls, which 
can be neglected as they are located "behind" the dominant walls. 

In the /9-space picture, the free motion before a collision is described by a 
null straight line. The effect of a collision on a particular wall wa(P) is easily 
obtained by solving (3.41), or (3.42), with only one term in the sum. This 
dynamics is exactly integrable: it suffices to decompose the motion of the 
/9-particle into two (linear) components: (i) the component parallel to the 
(timelike) wall hyperplane, and (ii) the orthogonal component. One easily 
finds that the parallel motion is left unperturbed by the presence of the wall, 
while the orthogonal motion suffers a (one-dimensional) reflection, with a 
change of the sign of the outgoing orthogonal velocity with respect to the 
ingoing one. The net effect of the collision on a certain wall w(/3) then is 
to change the ingoing velocity vector f M = df3^ jdr entering the ingoing free 
motion into an outgoing velocity vector v >iX given (in any linear frame) by 
the usual formula for a geometric reflection in the hyperplane w(/3) = 0: 

(w ■ w 

Here, all scalar products, and index raisings, are done with the /5-space metric 
G^y. Note that the collision law (5.3) leaves invariant the (Minkowski) length 
of the vector v^. Because the dominant walls are timelike, the geometric 
reflections that the velocities undergo during a collision, are elements of the 
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(5.3) 



orthochronous Lorentz group. Each reflection preserves the norm and the 
time-orientation; hence, the velocity vector remains null and future-oriented. 

From this perspective, we can also better understand the relevance of 
walls which are not timelike. Lightlike walls (like some of the subleading 
gravitational walls) can never cause reflections because in order to hit them 
the billiard ball would have to move at superluminal speeds in violation of the 
Hamiltonian constraint. The effect of spacelike walls (like the cosmological 
constant wall) is again different: they are either irrelevant (if they are "be- 
hind the motion"), or otherwise they reverse the time-orientation inducing a 
motion towards increasing spatial volume ("bounce"). 

The hyperbolic billiard is obtained from the /3-space picture by a radial 
projection onto the unit hyperboloid of the piecewise straight motion in the 
polywedge defined by the walls. The straight motion thereby becomes a 
geodesic motion on hyperbolic space. The "cushions" of the hyperbolic bil- 
liard table are the intersections of the hyperplanes (5.1) with the unit hy- 
perboloid, such that the billiard motion is constrained to be in the region 
defined as the intersection of the half-spaces wa(P) > with the unit hyper- 
boloid. As we already emphasized, not all walls are relevant since some of the 
inequalities wa(P) > are implied by others [26]. Only the dominant wall 
forms, in terms of which all the other wall forms can be expressed as linear 
combinations with non-negative coefficients, are relevant for determining the 
billiard. Usually, these are the minimal symmetry walls and some of the p- 
form walls. The billiard region, as a subset of hyperbolic space, is in general 
non- compact because the cushions meet at infinity (i.e. at a cusp); in terms 
of the original scale factor variables (3, this means that the corresponding 
hyperplanes intersect on the lightcone. It is important that, even when the 
billiard is non-compact, the hyperbolic region can have finite volume. 

Given the action (2.1) with definite spacetime dimension, menu of fields 
and dilaton couplings, one can determine the relevant wall forms and compute 
the billiards. For generic initial conditions, we have the following results, as 
to which of the models (2.1) exhibit oscillatory behaviour (finite volume 
billiard) or Kasner-like behaviour (infinite volume billiard) 

• Pure gravity billiards have finite volume for spacetime dimension D < 
10 and infinite volume for spacetime dimension D > 11 [34]. This 
can be understood in terms of the underlying Kac-Moody algebra AE^ 
[27]: as shown there, the system is chaotic precisely if the underlying 
indefinite Kac-Moody algebra is hyperbolic. 
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• The billiard of gravity coupled to a dilaton always has infinite volume, 
hence exhibits Kasner-like behavior [7, 2, 30]. 

• If gravity is coupled to p-forms (with p ^ and p < D — 2) with- 
out a dilaton the corresponding billiard has a finite volume [25] . The 
most prominent example in this class is D = 11 supergravity, whereas 
vacuum gravity in 11 dimensions is Kasner-like. The 3-form is crucial 
for closing the billiard. Similarly, the Einstein-Maxwell system in four 
(in fact any number of) dimensions has a finite-volume billiard (see 
[48, 59, 70] for a discussion of four-dimensional homogeneous models 
with Maxwell fields exhibiting oscillatory behaviour). 

• The volume of the mixed Einstein-dilaton-p-form system depends on 
the dilaton couplings. For a given spacetime dimension D and a given 
menu of p-forms there exists a subcritical open domain T> in the space of 
the dilaton couplings, i.e. an open neighbourhood of the origin \ p = 
such that: (i) when the dilaton couplings A p belong to V the general 
behaviour is Kasner-like, but (ii) when the \ p do not belong to V 
the behaviour is oscillatory [24, 30]. For all the superstring models, 
the dilaton couplings do not belong to the subcritical domain and the 
billiard has finite volume. Note, however, that the superstring dilaton 
couplings are precisely "critical", i.e. on the borderline between the 
subcritical and the overcritical domain. 

As a note of caution let us point out that some indicators of chaos must 
be used with care in general relativity, because of reparametrization invari- 
ance, and in particular redefinitions of the time coordinate; see [20, 43] for a 
discussion of the original Bianchi IX model. In this respect, we refer to [29] 
for a discussion of the link between the various time coordinates used in the 
analysis above: t, r and T. 

The hyperbolic billiard description of the (3+l)-dimensional homoge- 
neous Bianchi IX system was first worked out by Chitre [19] and Misner 
[64]. It was subsequently generalized to inhomogeneous metrics in [54, 45]. 
The extension to higher dimensions with perfect fluid sources was considered 
in [56] , without symmetry walls. Exterior p-form sources were investigated in 
[46, 47] for special classes of metric and p-form configurations. The uniform 
approach (based on the Iwasawa decomposition) summarized above comes 
from [29]. 
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6 Kac-Moody theoretic formulation: the E w 

case 



Although the billiard description holds for all systems governed by the ac- 
tion (2.1), the billiard in general has no notable regularity property. In 
particular, the dihedral angles between the faces, which can depend on the 
(continuous) dilaton couplings, need not be integer submultiples of ir. In 
some instances, however, the billiard can be identified with the fundamental 
Weyl chamber of a symmetrizable Kac-Moody (or KM) algebra of indefi- 
nite type 10 , with Lorentzian signature metric [26, 27, 23]. Such billiards are 
called "Kac-Moody billiards" . More specifically, in [26] , superstring models 
were considered and the rank 10 KM algebras E w and BE W were shown to 
emerge, in line with earlier conjectures made in [50, 51] n . This result was 
further extended to pure gravity in any number of spacetime dimensions, for 
which the relevant KM algebra is AE d , and it was understood that chaos 
(finite volume of the billiard) is equivalent to hyperbolicity of the underlying 
Kac-Moody algebra [27]. For pure gravity in D = 4 the relevant algebra 
is the hyperbolic algebra AE% first investigated in [38]. Further examples 
of emergence of Lorentzian Kac-Moody algebras, based on the models of 
[15, 21], are given in [23]. See also [32] and [40]. 

The main feature of the gravitational billiards that can be associated 
with KM algebras is that there exists a group theoretical interpretation of 
the billiard motion: the asymptotic BKL dynamics is equivalent (in a sense to 
be made precise below), at each spatial point, to the asymptotic dynamics of 
a one-dimensional nonlinear a-model based on a certain infinite dimensional 
coset space G/K, where the KM group G and its maximal compact subgroup 
K depend on the specific model. As we have seen, the walls that determine 
the billiards are the dominant walls. For KM billiards, they correspond 
to the simple roots of the KM algebra. As we discuss below, some of the 
subdominant walls also have an algebraic interpretation in terms of higher- 
height positive roots. This enables one to go beyond the BKL limit and to 
see the beginnings of a possible identification of the dynamics of the scale 
factors and of all the remaining variables with that of a non-linear a-model 

10 Throughout this section, we will use the abbreviations "KM" for "Kac-Moody", and 
"CSA" for Cartan subalgebra. 

11 Note that the Weyl groups of the E- family have been discussed in a similar vein in 
the context of [/-duality [61, 65, 3]. 
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defined on the cosets of the Kac-Moody group divided its maximal compact 
subgroup [28, 29]. 

The KM theoretic reformulation not only enables us to give a unified 
group theoretical derivation of the different types of walls discussed in the 
preceding section, but also shows that the (3 -space of logarithmic scale factors, 
in which the billard motion takes place, can be identified with the Cartan 
subalgebra of the underlying indefinite Kac-Moody algebra. The various types 
of walls can thus be understood directly as arising from the large field limit of 
the corresponding cr-models. It is the presence of gravity, which comes with 
a metric in scale-factor space of Lorentzian signature, which forces us to 
consider infinite dimensional groups if we want to recover all the walls found 
in our previous analysis, and this is the main reason we need the theory of 
KM algebras. For finite dimensional Lie algebras we obtain only a subset of 
the walls: one of the cushions of the associated billiard is missing, and one 
always ends up with a monotonic Kasner-type behavior in the limit t — > + . 
The absence of chaotic oscillations for models based on finite dimensional Lie 
groups is consistent with the classical integrability of these models. While 
they remain formally integrable for infinite dimensional KM groups, one can 
understand the chaotic behavior as resulting from the projection of a motion 
in an infinite dimensional space onto a finite dimensional subspace. 

For concreteness, we shall only consider one specific example here: the 
relation between the cosmo logical evolution of D = 11 supergravity and 
a null geodesic on E W /K(E W ) [28]. See Refs. [57, 58] for the extension 
of this correspondence to the two supergravities in D = 10, and Ref. [31] 
for a complementary check of the correspondence in the case of D = 11 
supergravity. We refer to [29] for a more general discussion of the link between 
KM billiards and KM coset models (including a discussion of the AE 3 case 
relevant for pure Einstein gravity in 3 + 1 dimensions). Note also that the 
relevance of non-linear cr-models for uncovering the symmetries of M-theory 
has also been discussed from a different, spacetime-covariant point of view 
in [71, 68, 69], but there it is E u rather than E w that has been proposed as 
a fundamental symmetry. 

The action defining the bosonic part of D = 11 supergravity reads 



S 
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d 
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where the spacetime indices a, (3, ... — 0,1, ... ,10, where £ 01 - 10 = -fl ; and 
where the four-form T is the exterior derivative of A: T = dA. Note the 
presence of the Chern-Simons term J-'AJ-'AA in the action (6.1). Introducing 
a zero-shift slicing (N l = 0) of the eleven-dimensional spacetime, and a time- 
independent spatial zehnbein 9 a (x) = E a i(x)dx\ the metric and four form 
T = dA become 

ds 2 = G Q/3 dx a dx p = -N 2 {dx ) 2 + G ab e a 6 b (6.2) 
T = ^F Oabc dx o A6 a A6 b A6 c + ^f abcd 6 a A6 b A6 c A0 d . 

We choose the time coordinate x° so that the lapse N = y/G, with G := 
detG a fe (note that x° is not the proper time 12 T = J Ndx°; rather, x° — > oo 
as T — > 0). In this frame the complete evolution equations of D = 11 
supergravity read 

ol^Cr OoLrcb) — gCry- J-b^s — -fr^J~ J- a /3-yS0 b 
-2GR a b (T,C) 

d( /~t /r-0a6c\ 1 r .abcaia-2a a bi 626364 -77 -77 

J — ^"0010203^61626364 

+ ^ G:f de[ab C c] de _ Q C e de jrdabc _ Q^Qjrdabc^ 
abed = 6JF ] e [abC e cd] + 49[ a ^ r 3 ft C( 2] (6.3) 

where a, 6 G {1, . . . , 10} and a, (3 G {0, 1, . . . , 10}, and i? a t(r, C) denotes 
the spatial Ricci tensor; the (frame) connection components are given by 
2G a dT d bc = Cabc+Cbca—C ca b+dbG ca +d c G a b—d a G bc withC^c = G ad Cd bc being 
the structure coefficients of the zehnbein d6 a = \C a bc Q h AQ c . [Note the change 
in sign convention here compared to above.] The frame derivative is d a = 
E l a (x)di (with E a iE l b = 5%). To determine the solution at any given spatial 
point x requires knowledge of an infinite tower of spatial gradients: one 
should thus augment (6.3) by evolution equations for d a Gb c , <9 a jF obcd , d a T bc de-, 
etc., which in turn would involve higher and higher spatial gradients. 

The geodesic Lagrangian on E w / K(E W ) is defined by generalizing the 
standard Lagrangian on a finite dimensional coset space G/K, where K is 
a maximal compact subgroup of the Lie group G. All the elements entering 
the construction of C have natural generalizations to the case where G is 
the group obtained by exponentiation of a hyperbolic KM algebra. We refer 

12 In this section we denote the proper time by T to keep the variable t for denoting the 
parameter of the one-dimensional er-model introduced below. 
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readers to [52] for basic definitions and results of the theory of KM algebras, 
and here only recall that a KM algebra g = q(A) is generally defined by 
means of a Cartan matrix A and a set of generators {e^, fi, hi} and relations 
(Chevalley-Serre presentation), where i,j = 1, . . . ,r = rankg(A). The ele- 
ments {hi} span the Cartan subalgebra (CSA) f), while the and fi generate 
an infinite tower of raising and lowering operators, respectively. The "maxi- 
mal compact" subalgebra t is defined as the subalgebra of q(A) left invariant 
under the Chevalley involution uu(hi) = —hi , a>(ej) = —fi , tv(fi) = — e*. 
In other words, t is spanned by the "antisymmetric" elements E a>s — E^ s , 
where E^ s = —u(E a ^ s ) is the "transpose" of some multiple commutator E atS 
of the ej's associated with the root a (i.e. [h,E a>s ] = a(h)E a ^ s for h £ h). 
Here s — 1, . . . mult(o;) labels the different elements of q(A) having the same 
root a. The (integer- valued) Cartan matrix of E\q is encoded in its Dynkin 
diagram. See top diagram in Fig. 2. 

The a-model is formulated in terms of a one-parameter dependent group 
element V = V(t) G E 10 and its Lie algebra valued derivative 

v{ t) '= ^7 V_1 W G e io = Lie ^io- (6-4) 
at 

The action is f dtC with 

C := n(t) _1 (v sym (t)\v sym (t)) (6.5) 

with a "lapse" function n(t) (not to be confused with N), whose variation 
gives rise to the Hamiltonian constraint ensuring that the trajectory is a 
null geodesic. The "symmetric" projection t> sym := \{v + v T ) eliminates the 
component of v corresponding to a displacement "along fi" , thereby defining 
an evolution on the coset space E w / K(E W ). (.|.) is the standard invariant 
bilinear form on the KM algebra [52]. 

Because no closed form construction exists for the raising operators E a>s , 
nor their invariant scalar products {E a>s \E^ = N^ t 5^ +/3 , a recursive ap- 
proach based on the decomposition of Eiq into irreducible representations of 
its SL(10, M) subgroup was devised in [28]. [See also [57] for a decomposition 
of E w based on its D Q subgroup, and [58] for a decomposition of E 10 based on 
its A 8 x Ai subgroup.] Let a±, . . . , a 9 be the nine simple roots of A 9 = sZ(10) 
corresponding to the horizontal line in the E w Dynkin diagram, and a the 
"exceptional" root connected to a 3 . [This root is labelled 10 in Fig.2.] Its 
dual CSA element h enlarges A 9 to the Lie algebra of GL(10). Any positive 
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root of E w can be written as 



9 



a = £a + mPotj (£, m? > 0). 



(6.6) 



We call i = i(a) the "level" of the root a. This definition differs from the 
usual one, where the (affine) level is identified with m 9 and thus counts the 
number of appearances of the over-extended root ag in a [38, 53]. Hence, our 
decomposition corresponds to a slicing (or "grading") of the forward light- 
cone in the root lattice by spacelike hyperplanes, with only finitely many 
roots in each slice, as opposed to the lightlike slicing for the Eg representa- 
tions (involving not only infinitely many roots but also infinitely many affine 
representations for m 9 > 2 [38, 53]). 

The adjoint action of the Ag subalgebra leaves the level £(a) invariant. 
The set of generators corresponding to a given level i can therefore be de- 
composed into a (finite) number of irreducible representations of Ag. The 
multiplicity of a as a root of Eiq is thus equal to the sum of its multiplicities 
as a weight occurring in the SL(10, M.) representations. Each irreducible rep- 
resentation of Ag can be characterized by its highest weight A, or equivalently 
by its Dynkin labels (pi, . . . ,p 9 ) where Pfc(A) := («&, A) > is the number 
of columns with k boxes in the associated Young tableau. For instance, the 
Dynkin labels (001000000) correspond to a Young tableau consisting of one 
column with three boxes, i.e. the antisymmetric tensor with three indices. 
The Dynkin labels are related to the 9-tuple of integers (m 1 , . . . , m 9 ) appear- 
ing in (6.6) (for the highest weight A = —a) by 



where is the inverse Cartan matrix of Ag. This relation strongly constrains 
the representations that can appear at level £, because the entries of are 
all positive, and the 9-tuples (pi, . . . ,pg) and (mi, . . . , m 9 ) must both consist 
of non-negative integers. In addition to satisfying the Diophantine equations 
(6.7), the highest weights must be roots of E w , which implies the inequality 
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(6.7) 



9 




(6.8) 
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All representations occurring at level £ + 1 are contained in the product of 
the level-£ representations with the £ = 1 representation. Imposing the Dio- 
phantine inequalities (6.7), (6.8) allows one to discard many representations 
appearing in this product. The problem of finding a completely explicit and 
manageable representation of E w in terms of an infinite tower of Ag represen- 
tations is thereby reduced to the problem of determining the outer multiplic- 
ities of the surviving Ag representations, namely the number of times each 
representation appears at a given level I. The Dynkin labels (all appearing 
with outer multiplicity one) for the first six levels of E\q are 



i= 1 


(001000000) 




1 = 2 


(000001000) 




£ = 3 


(100000010) 




£ = A 


(001000001) , 


(200000000) 


£ = 5 


(000001001) , 


(100100000) 


£ = 6 


(100000011) , 


(010001000) , 




(100000100) , 


(000000010) . 



(6.9) 

The level £ < 4 representations can be easily determined by comparison 
with the decomposition of E$ under its Aj subalgebra and use of the Jacobi 
identity, which eliminates the representations (000000001) at level three and 
(010000000) at level four. By use of a computer and the E 10 root multiplici- 
ties listed in [53, 4], the calculation can be carried much further [39]. 

From (6.9) we can now directly read off the GL(10) tensors making up 
the low level elements of Eiq. At level zero, we have the GL(10) gener- 
ators K a b obeying [K a b , K c d ] = K a d 5 b — K c b 5%. The eio elements at lev- 
els £ = 1,2,3 are the GL(10) tensors E aia * a3 , E ai - a(i and E a °\ ai - a * with 
the symmetries implied by the Dynkin labels. The cr-model associates to 
these generators a corresponding tower of "fields" (depending only on the 
"time" t): a zehnbein h a b (t) at level zero, a three form A abc (t) at level 
one, a six-form A ai _,_ a6 (t) at level two, a Young-tensor A ao \ ai ... as (t) at level 3, 
etc. Writing the generic E w group element in Borel (triangular) gauge as 
V(f) = expX h (t) -expX A (t) with X h {t) = h\K b a and X A (t) = ±A abc E abc + 
±A ai .„ a6 E ai - a6 + ^A ao \ ai ... a& E a ^- a * + ..., and using the E w commutation 
relations in GL(10) form together with the bilinear form for Eiq, we find up 
to third order in level 

nC = \(g<"g M -g<*g<*)grtgcd+Ul D ^i™ DAaiaaaa 
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+^DA ai ... a6 DA a ^ + \±DA Mai ... a ,DA a ^--- a » (6.10) 

where g ab = e a c e b c with e a b = (exp h) a b , and all "contravariant indices" have 
been raised by g ab . The "covariant" time derivatives are defined by (with 
dA = A) 

DA 



= dA aia2a3 (6.11) 
dA ai a6 + \0A^ aia2a:i dA a4a5a6 ^ 
dA ai |d 2 ...ag ~t" 4 :< ^A^ aia2a-3 0A ai a9 ^ 
A1dA^ aia2a . A A ai ... a9 ^ -|- 280^4^ aia2a3 j4 a4asa6 c}j4 a7aga g^ . 



DA ai ^ a& 

DA ai |a 2 ... ag 



Here antisymmetrization [...], and projection on the £ = 3 representation 
(...), are normalized with strength one (e.g. [[...]] = [•••])■ Modulo field 
redefinitions, all numerical coefficients in (6.10) and (6.11) are uniquely fixed 
by the structure of Eiq. See also [31] for a different perspective on the 
equations of motion of the E 10 / K(E 10 ) a-model. 

The Lagrangian (6.5) is invariant under a nonlinear realization of E 10 
such that V(t) — > k g (t)V(t)g with g e E 10 ; the compensating "rotation" 
k g (t) being, in general, required to restore the "triangular gauge". When 
g belongs to the nilpotent subgroup generated by the E abc , etc., this sym- 
metry reduces to the rather obvious "shift" symmetries of (6.10) and no 
compensating rotation is needed. The latter are, however, required for the 
transformations generated by F abc = (E abc ) T , etc. The associated infinite 
number of conserved (Noether) charges are formally given by J = Ai~ 1 dAi, 
where M. = V T V. This can be formally solved in closed form as 

M(t)=M(0)-esxp(tJ). (6.12) 

The compatibility between (6.12) (indicative of the integrability of (6.10)) 
and the chaotic behavior of g a b(t) near a spacelike singularity is discussed in 
[29]. 

The main result of concern here is the following: there exists a map 
between geometrical quantities constructed at a given spatial point x from 
the supergravity fields G^(x°,a;) and A^ up (x°,x) and the one-parameter- 
dependent quantities g a b(t), A a b c (t), . . . entering the coset Lagrangian (6.10), 
under which the supergravity equations of motion (6.3) become equivalent, 
up to 30th order in height, to the Euler-Lagrange equations of (6.10). In the 
gauge (6.2) this map is defined by t = x° = J dT/^/G and 

9ab{t) = G ab (t,x) (6.13) 
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-DA aia2a3 (t) J- 3 aifl2 a3 (t, X) 

DA a ^(t) = -^ ai - a6blfe2b364 ^ 1&2& 3 64 (t,x) 

The expansion in height ht(a) = £ + ^m J , which controls the iterative 
validity of this equivalence, is as follows: the Hamiltonian constraint of 
the coset model (6.10) contains an infinite series of exponential coefficients 
exp ( — 2a(/3)), where a runs over all positive roots of Eiq, and where 
[3 a = — h a a parametrize the CSA of E w . The billiard picture discussed above 
is equivalent to saying that, near a spacelike singularity (t — > oo), the dy- 
namics of the supergravity fields and of truncated versions of the E 10 coset 
fields is asymptotically dominated by the (hyperbolic) Toda model defined 
by keeping only the exponentials involving the simple roots of Eiq. Higher 
roots introduce smaller and smaller corrections as t increases. The "height 
expansion" of the equations of motion is then technically defined as a for- 
mal BKL-like expansion that corresponds to such an expansion in decreasing 
exponentials of the Hamiltonian constraint. On the supergravity side, this 
expansion amounts to an expansion in gradients of the fields in appropriate 
frames. Level one corresponds to the simplest one-dimensional reduction of 
(6.3), obtained by assuming that both and Ax^ depend only on time 
[25]; levels 2 and 3 correspond to configurations of G^ u and A\^ v with a 
more general, but still very restricted x-dependence, so that e.g. the frame 
derivatives of the electromagnetic field in (6.3) drop out [33]. In [28, 31] it is 
checked that, when neglecting terms corresponding to ht(a) > 30, the map 
(6.13) provides a perfect match between the supergravity evolution equations 
(6.3) and the E w coset ones, as well as between the associated Hamiltonian 
constraints. (In fact, the matching extends to all real roots of level < 3.) Let 
us also mention in passing (from [31]) that the E w coset action is not compat- 
ible with the addition of an eleven-dimensional cosmological constant in the 
supergravity action (an addition which has been proven to be incompatible 
with supersymmetry in [5]). 

It is natural to view the map (6.13) as embedded in a hierarchical sequence 
of maps involving more and more spatial gradients of the basic supergravity 
fields. The height expansion would then be a way of revealing step by step 
a hidden hyperbolic symmetry, implying the existence of a huge non-local 
symmetry of Einstein's theory and its generalizations. Although the validity 
of this conjecture remains to be established, one can at least show that there 
is "enough room" in E 10 for all the spatial gradients. Namely, the search 
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for affine roots (with m 9 = 0) in (6.7) and (6.8) reveals three infinite sets 
of admissible A 9 Dynkin labels (OOlOOOOOn), (OOOOOlOOn) and (lOOOOOOln) 
with highest weights obeying A 2 = 2, at levels t — 3n + 1, 3n + 2 and 3n + 3, 
respectively. These correspond to three infinite towers of e elements 

E a , a blb2b3 , E a , a bl - b6 , E a , a bolbl - bs (6.14) 

u, \ . ..(-in ' (*l . . .tin ' u,\. ..u, n \ / 

which are symmetric in the lower indices and all appear with outer multi- 
plicity one (together with three transposed towers). Restricting the indices 
to cij = 1 and bi G {2, 10} and using the decomposition 248—^80+84+84 
of E% under its SL(9) subgroup one easily recovers the affine subalgebra 
Eg C E w . The appearance of higher order dual potentials (d la Geroch) in 
the .Eg-based linear system for D = 2 supergravity [14] indeed suggests that 
we associate the E w Lie algebra elements (6.14) to the higher order spa- 
tial gradients d ai ■ ■ ■ d an A blb2ba ,d ai ■ ■ ■ d an A bl ^ bfj and d ai ■ ■ ■ d an A bo \ bl _ bs or to 
some of their non-local equivalents. Finally, we refer to [29] for a more ge- 
neral discussion of the height expansion of Kac-Moody cr-models, and we note 
that the approach of [28, 29] can be extended to other physically relevant 
KM algebras, such as BE W [26, 21] and AE n [27]. 
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